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Luttinger liquids with multiple Fermi edges: Generalized Fisher-Hartwig conjecture 
and numerical analysis of Toeplitz determinants 
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It has been shown that solutions of a number of many-body problems out of equilibrium can be 
expressed in terms of Toeplitz determinants with Fisher-Hartwig (FH) singularities. In the present 
paper, such Toeplitz determinants are studied numerically. Results of our numerical calculations 
fully agree with the FH conjecture in an extended form that includes a summation over all FH 
representations (corresponding to different branches of the logarithms). As specific applications, 
we consider problems of Fermi edge singularity and tunneling spectroscopy of Luttinger liquid with 
multiple-step energy distribution functions, including the case of population inversion. In the energy 
representation, a sum over FH branches produces power-law singularities at multiple edges. 
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I. INTRODUCTION 

For more than half a century, quantum many-body sys- 
tems remain one of central research directions in the con- 
densed matter physics. There is a number of quantum 
many-body problems that are of fundamental physical 
importance and, at the same time, possess an exact solu- 
tion. These are the Anderson orthogonality catastrophe^, 
Fermi edge singularity^ (FES), Luttinger liquid^ (LL) 
zero-bias anomaly^-, and Kondo problems^. It has been 
realized long ago that these problems are, in fact, deeply 
interconnected, both from the point of view of the un- 
derlying physics and of the mathematics involved. Such 
connections have been used, e.g., for the representation 
of the dynamics of the Kondo problem as an infinite 
sequence of Fermi-edge-singularity events^. These rela- 
tionships betvifeen many-body problems extends beoynd 
fermions and encompass also interacting bosons (e.g., the 
Lieb-Liniger modeli), one-dimensional Heisenberg chains 
and etciS 

In recent works by two of us with Gefen2r2i, non- 
equilibrium realizations of some of these problems have 
been investigated. For this purpose, we have devel- 
oped a non-equilibrium bosonization technique general- 
izing the conventional bosonizationi2r2& onto problems 
with non-equilibrium distribution functions. We have 
shown that the relevant correlation functions can be ex- 
pressed through Fredholm determinants of "counting" 
operator. The information on the specific type of the 
problem, as well as on different aspects of the interac- 
tion, is encoded in the time-dependent scattering phase 
of the counting operator. The findings of Refs. MWu have 
demonstrated that the above classical many-body prob- 
lems are even more closely connected than has been pre- 
viously understood, extending the interrelations into the 
non-equilibrium regime. 

The "counting" operators governing the simplest (one- 



particle) non-equilibrium Green functions in the above 
models can be reduced to Toeplitz matrix form upon 
regularization and discretization^^. The electron energy 
distribution function then determines the symbol of the 
Toeplitz matrix. The most interesting situation arises 
when the distribution function has multiple steps ( "Fermi 
edges"), which results in step- like singularities of the 
symbol. According to the Fisher-Hartwig conjectureii, 
this leads to a non-trivial power-law behavior of the 
correlation functions. Recent progress in the analysis 
of Toeplitz determinants with Fisher-Hartwig singular- 
ities has allowed to establish their leading asymptotic 
behavior!^. In Ref. |Tl| a generalized Fisher-Hartwig 
conjecture was put forward that includes a summation 
over all Fisher-Hartwig representations (corresponding 
to different branches of the logarithm of the symbol). 
This yields also terms with subleading power-law factors. 
While these terms are formally smaller (as compared to 
the leading term) when one considers the Green function 
in the time representation, they contain different oscilla- 
tory exponents. Therefore, after a transformation to the 
energy representation, they produce power-law singular- 
ities at different edges, which makes these terms phys- 
ically important. The extended version of the Fisher- 
Hartwig conjecture is also expected to be of interest from 
the purely mathematical point of view. 

In the present paper we perform a numerical analysis 
of Toeplitz determinants with Fisher-Hartwig singulari- 
ties. The results of numerical calculation fully confirm 
the extended conjecture for the asymptotic (long-time or 
low-energy) behavior. Furthermore, the numerics allows 
us to explore correlation functions in the entire energy 
range. To be specific, we focus on two fermionic prob- 
lems: (i) the Fermi-edge singularity in X-ray absorption 
and (ii) the tunneling density of states (TDOS) of a non- 
equilibrium Luttinger liquid. 

The structure of the present paper is as follows. Sec- 



tion|TT] contains a brief review of the connection between 
one-particle correlation functions of many-body problems 
and Toeplitz determinants. In Sec. lIIII we present the ex- 
tended version of the Fisher-Hartwig conjecture, as well 
as illustrate it and discuss its implucations on examples 
relevant to our many-body problems. In Sec. lIVI we calcu- 
late the Toeplitz determinants (and thus the correlation 
functions under interest) numerically and compare the 
exact results with the asymptotic formulas. Our find- 
ings are summarized in Section |Vl where we also discuss 
prospects for future research. 
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II. MANY-PARTICLE PROBLEMS AS 
FREDHOLM DETERMINANTS 



FIG. 1: Rectangular shaped pulse in the phase S{t). 



A. Fermi edge singularity 



B. Luttinger liquid: tunneling spectroscopy 



The FES problem describes the scattering of conduc- 
tion electrons off a localized hole which is left behind by 
an electron excited into the conduction band. Histori- 
cally, the FES problem was first solved by exact sum- 
mation of an infinite diagrammatic series^. While in the 
FES problem there is no interaction between electrons in 
the conducting band, it has many features characteris- 
tic of genuine many-body physics. Despite the fact that 
conventional experimental realizations of FES are three- 
dimensional, the problem can be reduced (due to the local 
and isotropic character of the interaction with the core 
hole) to that of one-dimensional chiral fermions. For this 
reason, bosonization technique can be effectively applied, 
leading to an alternative and very elegant solution^^. 

One can consider the FES out of equilibrium^^, with 
an arbitrary electron distribution function n{e). This 
problem can be solved within the framework of non- 
equilibrium bosonization^, with the following results for 
the emission/absorption rates: 
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AA^(27r-2i5o) 
2ttv{1 ± iAr)(i-'5"A)' 



(1) 



Here Sqis the s-wave electronic phase shift due to the scat- 
tering of conduction electrons off the core hole. Further, 
At- [27r — 26o\ is the Fredholm determinant (normalized to 
its value at zero temperature) 
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AAS] 



dct[l + (e 



-iS 



l)h] 



Ar[S,T = 0] dct[l + (e-^-^ - l)h{T = 0)] 



(2) 

The phase 6 is an operator local in time t conjugate to 
electron energy e and has characteristic rectangular shape 

(Fig. m 



S{t) = 6 [9{t) - 9{t - t)] . 



(3) 



The connection of the non-equilibrium FES problem to 
Fredholm determinants is summarized in the first row of 
Table H 



The tunneling spectroscopy technique allows one to ex- 
plore experimentally Keldysh Green functions of an inter- 
acting system that carry information about both tunnel- 
ing density of states and energy distribution. Recent ex- 
periments on carbon nanotubes and quantum Hall edges 
have proved the efficiency of this technique in the con- 
text of ID systems2ii2^. The technological and experi- 
mental advances motivate the theoretical interest in the 
tunneling spectroscopy of strongly correlated ID struc- 
tures away from equilibriunt^i^^^— . 

In the case of a LL formed by ID interacting fermions, 
the Keldysh Green function may be evaluated theoreti- 
cally via the non-equilibrium bosonization technique. As- 
suming that a long LL conductor is adiabatically cou- 
pled to two reservoirs (modeled as non-interacting ID 
wirea^"— ) with distribution functions nfl(e) and ^^(e) 
respectively, one obtains for the Green functions of the 
right movers^ 
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27ru (l±iAr)i+T ' 

where u = v/K is the sound velocity, 

j=il-K)y2K, (5) 

and 

K ^(l+g/TTv)-^^^ (6) 

is the standard LL parameter in the interacting region. 
The determinants A^7-[(5^] (77 = R, L) are given by Eq. 
([2]) with n{e) replaced by the corresponding distribution 
functions n^(e) and 



S = Sri = TT 



l + f]K 
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(7) 



The connection of the Luttinger liquid Green functions 
to Fredholm determinants is summarized in the last two 
rows of Table ID 





5r 


5l 


7 


G^^sir) 


2(7r - 5o) 





5g 25n 

7r2 TT 


Glir) 


2^/K 


2"^ 


2K 


Glir) 




^-ijf 


2K 



TABLE I: Non-equilibrium correlation functions of many- 
body problems: Fermi edge singularity (Gpgg(r)), Green 
functions of right- and left-moving fermions in a LL {G^{t) 
and Gl (t)). Evaluation of these correlation functions yields 
the results in the form of Fredholm-Toeplitz determinants 
T~'^~^ Ar[5r]Al[5l]. The corresponding phases Sr^l are pre- 
sented in the second and third columns. (For LL an adiabatic 
coupling to reservoirs on the scale of the characteristic plas- 
mon wave length is assumed.) The determinants are normal- 
ized to their values at zero temperature. The exponent 7 
governing the zero-temperature correlation function is shown 
in the last column. 



It is worth emphasizing that the the rectangular shape 
Q of the pulse with the amplitude ([7]) is valid in the 
case when the coupling to reservoirs is smooth on the 
scale of the plasmon wave length v/T, u/T. In the op- 
posite regime the pulse S{t) entering Q is fractionalized 
in a sequence of rectangular pulses^. In the long-wire 
limit the corresponding determinant splits into a prod- 
uct of single-pulse (i.e Toeplitz-type) determinants. For 
definiteness, we focus on the adiabatic case in this paper. 



C. Ultraviolet regularization and reduction to 
Toeplitz matrix 

Due to characteristic rectangular shape (|31) of the 
pulses 6{t) the Fredholm determinants At-((5) are in fact 
of the Toeplitz form. Specifically, one can write 



A^[^] = det[l + P{e-'^ - l)hP] . 
Here we have defined the projection operator 
y{t), for te [0,t] 

, otherwise . 
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(8) 
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The form ([8]) is convenient for peforming the ultraviolet 
regularization of the determinant A,- [S] . Specifically, we 
discretize the time t by introducing an elementary time 
step At = tt/A, such that tj — jAt. This corresponds 
to restricting the energy variable e to the range [—A, A]. 
We arrive then at a finite-dimensional determinant 



An[S] - det[/(tj - tk)] , < i, fc < iV - 1 . 



(10) 



Here N = tA/tt and f{tj — tk) is Fourier transform of 
the function 



The matrix elements f{tj — tk) depend on j and k via 
the difference j — k only, so that the obtained matrix is 
of Toeplitz type. 

In order to bring Eqs. (ITU)) . ([TT]) to the canonical form 
used in the theory of Toeplitz matrices, we have to de- 
fine the function /(e) on the unit circle \z\ — 1. This 
is easily done by identifying the polar angle 6 £ [— tt, tt] 
parametrizing the unit circle via z = e'^ with the appro- 
priately rescaled energy: 9 = ire/ A. However, if this is 
done directly with the function (jllI) . a non-physical jump 
will arise ed, = ±tt. In order to eliminate it, one has to 
introduce an additional phase factor into the definition 
of /(e): 



/(e)-[l + n(e)(e- 



-iS 



l)]e- 



(12) 



/(e) = l + n(e)(e- 
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(11) 



After the mapping to the unit circle, z — e''^^/'^, this de- 
fines a function f{z) (known as the symbol of the Toeplitz 
matrix) that is perfectly smooth at z = — 1. It will, 
however, have discontinuities ( "Fisher-Hartwig singular- 
ities") at the positions z = e*'^'^^/^ if the distribution 
function n(e) has such discontinuities ("Fermi edges") at 
e j . We will be interested in the situation when there are 
several (at least two) such discontinuities. 

It is worth emphasizing that the regularization (jlOp . 
(1121) makes explicit the dependence of the determinant 
At- ((5) on the integer part of 5/2tt (thus making re- 
dundant the procedure of analytical continuation from 
6 G [— 7r,7r] to larger \S\ discussed in Ref. 9). This al- 
lows us to directly compute the determinant at arbitrar- 
ily large (by absolute value) 5. 

As the matrix {f{tj — tk)} with < j, fc < A^ — 1 is 
of Toeplitz form, results concerning the large- A'^ asymp- 
totic behavior of Toeplitz determinats Ajv can be ap- 
plied. We summarize them in the next section. Physi- 
cally, the large- A^ limit corresponds to the regime of long 
time T, i.e. to infrared asymptotics of correlation func- 
tions under interest. (In the energy representation, this 
translates into low-energy behavior around singularities.) 

Furthermore, Eqs. (|10p . ([T^ are also very convenient 
for numerical evaluation of the determinant Ar{d), pro- 
viding us access to the full time (or, after Fourier trans- 
formation, energy) dependence of the correlation func- 
tions. 



III. ASYMPTOTIC PROPERTIES OF 
TOEPLITZ DETERMINANTS 

Toeplitz matrices and operators were introduced by 
O. Toeplitz a century ago. Since this time, asymptotic 
properties of Toeplitz determinants have been in a fo- 
cus of interest of mathematicians, starting from the 1915 
paper'^^ that was the first research paper by G. Szego. 
The Szego theorems^ valid for a smooth symbol yield 
the large- Af asymptotics of the determinant, which is ex- 
ponential in A^, with an A^-independent prefactor. As was 
realized by Fisher and Hartwigii, in the case of a symbol 
with singularities, the asymptotics acquires, in addition 



to the exponential factor, also a power-law factor. Thus, 
the infrared behavior of the Toeplitz determinant pUj) in- 
cludes non-trivial power-law factors if the function f{z) 
is not smooth on the unit circle. The simplest example 
is the zero-temperature determinant^^ 
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that has a power-law dependence on time in the long-time 
limit (At > 1). 



Let us now consider a distribution function with mul- 



tiple steps: (cf. Fig. [2]): 
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(14) 



where < a^ < 1, for j — 0, . . . , m. We are inter- 
ested in the Toeplitz determinant (|10|) for the multi- 
step distribution function (fT4| . Let us split the phase 
(5 as (5 = 2nM + 6', where M is integer and \6'\ < tt. 
We find it convenient to normalize the determinat by 
its zero-temperature value (|13p: the normalized determi- 
nant will be denoted as At-{S). According to the ex- 
tended Fisher-Hartwig conjectureii, it has the following 
long time asymptotics 



Ar{S) = 



exp j-iTfi' - t/2t^] 
G{l-5/2n)G{l + d/2n) 



na + ...+n„ 
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(1 + ...). 



(15) 
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Here we use the following notations: the exponents /3' 
(satisfying |Re/3^| < 1/2) are 

I3'i = — [ln(l - Qi+i + Qi+ie^^^) - ln(l - ai + a^e"*'')] , 

(16) 
the dephasing rate reads 






(17) 



G{x) is the Barnes G-function, Ujk = \ej — (k\, and 
^' = —^eJ2jPj^j- Note that the ultraviolet regular- 
ization A does not enter the normalized determinant. 
The asymptotic p^ is valid provided that rUjk ^ 1 
for all j y^ k. The summation goes over all sets of integer 
no, ... , rim satisfying no -t- ... -I- n,„ = — M; each such set 
yields the corresponding oscillatory exponent e*^ ^j "^"^^ . 
Equation (llSp presents explicitly the leading asymptotic 
behavior for the factor multiplying each of these expo- 
nents. Apart from this dominant term, there will be in 
general also subleading (in powers of 1/i) terms corre- 
sponding to the same exponent; these are abbreviated 
by -I- ... in the last bracket. 

The asymptotics (llSp has a long history. The form 
of its leading term (the one with the slowest decay 
in r, i.e. with the smallest exponent a (no , . . . n„i) — 



2 Re 



J2i>k PjPk ) was suggested back in 1968 by Fisher 



and Hartwig^^. Since then, significant efforts were in- 
vested into the exact formulation, the proof, and exten- 
sions of the Fisher-Hartwig conjecture. For the case when 
a unique combination of integers ni exists that minimizes 
the exponent a (no , . . . n^), the leading asymptotic term 
(including the corresponding numerical coefficient indi- 
cated in (lisp ) was rigorously derived by Ehrhardli^. In a 
recent seminal paper^® the theorem due to Ehrhardt was 
generalized for the case when there are several distinct 
sets of integers {rii] sharing the same minimal value of 
the exponent a (no , . . . Um)- It was proven that the lead- 
ing term of the asymptotic expansion of the determinant 
At- [(5] at large t is given by p3|) where the sum should 
be restricted to the sets {n^} minimizing the exponent 
a(no , . . .n„i). 

More recently, two of us and Gcfen^^ formulated and 
extended version of the Fisher-Hartwig conjecture [as 
shown in Eq. p^ ] that includes a sum over all sets 
{rii} (which correspond to different branches of the log- 
arithms) and captures the leading term of the expansion 
at every oscillation frequency '^^rijej. This extension 
is very natural from the point of view of continutiy, as, 
under change of parameters, the dominant branch (that 
determines the leading asymptotics given by Ref.llSl) may 
become subdominant. This is particularly transparent in 
the energy representation of our problem discussed be- 
low: different branches then correspond to singularities 
near different energies; clearly, such a singularity will per- 
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FIG. 2: Left panel. Triple-step distribution: an example of multi-step distribution (|14|) . 
distribution of the type (|14p with identically zero dephasing. 
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Simplest non-trivial 



sist even when its exponent vifill become subdominant 
with respect to a singularity at other energy. Further- 
more, the summation over branches has a clear physical 
meaning: it corresponds to processes including transfer of 
one or several electrons between different Fermi edges^i. 
To illustrate how Eq. P3|) works, let us consider a 
simple case of the determinant at the phase 6 = 47r, 
which can be evaluated exactly by a "refermionization" 
procedure^. We assume for definiteness that the distri- 
bution function is of double step form with jumps at eo 
and ei. (Generalization to a distribution function with 
more than two jumps is straightforward.) The exact re- 
sult reads 

A,[47r] = ai(ai - 1) [(ei - eo)^^ - 2] e-^(^o+.i)r 

+ (l-ai)2e-2"''" + a?e-2"i-. (ig) 



On the other hand, considering the expansion (|15|) one 
gets at (5 — >■ 47r 



00 = (ai - 1) 
I3'i = -ai 



27r 
2^ 



(19) 
(20) 



with d' — > 0. Observing now that G{x) has k-th order 
zero ai X = —k + I for any positive integer k we con- 
clude that at (5 = 47r (or generally for any S being integer 
multiple of 2tt) the sum in ((T5|) becomes finite. In the 
present case only the terms with (no, rii) — (—1,-1), 
(—2,0) and (0,-2) contribute, yielding 



A,(4^) ~ ai(ai - l)(ei - eo)V2e-(^«+^i)- 



+ {l-ai)^e- 



+ a-,e 



-2ieiT 



Comparing this asymptotic formula to the exact result 
(fT8)) . we see that Eq. ([2T|) indeed perfectly reproduce 
leading factors for each oscillation frequency. The only 
term missing in Eq. (j2ip is 



which represents a small correction (due to an additional 
factor ex r^^) to the leading term at the same frequency 
eo + ei, 



ai(ai-l)(ei-eo)V2e-^(^«+^i)^ 



(23) 



Such terms representing small power-law corrections to 
the leading contribution at the same frequency are indi- 
cated in Eq. P3|) by the symbol -(- — 

Let us note that, while being small with respect to the 
leading term at the same frequency, these corrections are 
not necessarily small with respect to leading terms at 
other frequencies. In particular, in the considered exam- 
ple the correction term on the frequency eo + ei, Eq. ([2^ 
is of the same order as the terms oscillating with frequen- 
cies 2eo and 2ei that are taken into account by Eq. (|21[) . 

Thus, Eq. (fTS]) captures explicitly the leading term 
for each frequency. A mathematically rigorous proof 
of this generalized form of the Fisher-Hartwig conjec- 
ture remains to be developed. Also, one may hope that 
it is possible to generalize Eg. dTSI) further, accounting 
also for sub-leading contribution (indicated as + . . . in 
Eq. ([TSl) '). A construction of such a full asymptotic ex- 
pansion of the Toeplitz determinant was discussed very 
recently in Ref. |3a for the special case /(e) = 1 + (e^** — 
1)0 {U — |e|), where Q{x) is the Heaviside theta function. 

It is worth mentioning that for 5 = 2ti Eq. ([Tsl) repro- 
duces the exact result 

A^(27r) = (1 - a)e-"°^ -f ae'''^^ (24) 

without any corrections at all. While Eq. (IM)) is written 
for a double-step distribution, this statement is valid for 
any multi-step distribution as well. The only non-zero 
terms in Eq. ()15p for ^ = 27r are those with all rij being 
equal to zero except for one equal to —1. The deter- 
(21) minant demonstrates oscillations at frequencies e^. All 



2ai(ai - l)e-^('^°+^l)^ 



(22) 



corrections of the type -I- ... in Eq. (|15p vanish. This im- 
plies that for values of (5 ~ 27r the correction terms + . . . 
in Eq. ([T5|) have additional smallness. 

Having clarified the status of the expansion ([15]) , let us 
now discuss its implications. In a generic case, the power- 
law decay of A,- (5) is cut off by the non-equilibrium de- 
phasing time Tfj, given by Eq. (J17p . Quite remarkably. 



the dephasing time is an oscillating function of the phase 
d which translates, e.g., into the non-monotonous de- 
pendence of T^ on the interaction strength in Luttinger 
liquid^. The dephasing is absent when 5' — in which 
case A^ (S) can be represented in terms of a free fermionic 
theory. 

Dephasing is also absent for the case when all Qj = 0, 1. 
This corresponds to the case of a pure electronic state 
(i.e. characterized by a wave function rather than by a 
density matrix). The simplest non-trivial distribution of 
the type Eq. p^ that has this property is the triple-step 
distribution of Fig. [5] (right panel). We stress that in 
this "ideal inverse population" case the dephasing rate 



is identically equal to zero, regardless of the value of the 
phase 6. Apart from being interesting on a pure the- 
oretical grounds, the distributions realizing the inverse 
population of electronic states are also expected to be ex- 
perimentally relevant, as they are inevitably generated in 
course of evolution of a smooth perturbation of electronic 
density if the spectral curvature is taken into accounlj^. 

The power-law decay of At-(5) in the time domain is 
translated into the singular energy dependence of cor- 
relation functions in energy representation. Specifically, 
every term in the expansion (|15p gives rise to a singular 
contribution 
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(25) 



to the Fourier transform of Ar[S]/{iT)''^^ with the expo- 
nent 



j-2j2PjPk 
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In Sec. II VI we will compare Eqs. (P5|). (|26p with the results 
of numerical evaluation of Toeplitz determinants. 



IV. NUMERICAL ANALYSIS 

In this Section, we present results of the numerical 
analysis of the Toeplitz determinants (|10p which allows us 
to evaluate the many-body Green functions in the whole 
range of times (energies). We will further demonstrate 
that the numerics gives full support to the asymptotic 
expansion ([25]), (P5| . 



A. Numerical procedure 

To be specific, we will consider fermions with the fol- 
lowing two types of many-step distributions: (i) double- 
step distribution 



1 = ao, 
nd{e) ^ { ai = 1/3, 
= a2 , 



e < eg = -C//3 

eo < e < ei = 2t//3 (27) 

ei < e, 



and (ii) triple-step distribution with the 
verse population (Fig. [21 right panel) 
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e < eo = -3C//4 
eo < e < £1 == -C//2 
ei < e < £2 = U/A 

£2 < £• 



(28) 



In these equations we have expressed all the energies e^ 
in terms of characteristic scale U — e„i — eo associated 
with the distribution function. 

Let us consider the normalized determinant A,- [[/, S] 
and its finite-dimensional approximation A]s[[U/A,S]. 
Here we made explicit the dependence of the determi- 
nants on U. At T and U fixed, An[U/A,S] has a finite 
limit as A — > oo which is a cutoff-independent function 
of the dimcnsionless variable Ut only: 



Ar[U.,S] = A[Ut,5] = lim A 



A— >oo 
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(29) 
(30) 



Equation (1301) constitutes the starting point for our 
numerical analysis. With a simple Mathematica code we 
are able to go within a quite short computation time 
up to the size of the Toeplitz matrix N = 500, which 
is typically sufficient for the convergence to the large- iV 
limit for relevant values of Ut. 

The convergence properties of our procedure become 
generally worse at large 6. Thus, we chose to illustrate 
them with the calculation of the determinant at the phase 
47r — 0.6 which is larger then any phase we will encounter 
in the next section. This choice also enables us to demon- 
strate clearly the presence of the correction terms indi- 
cated by -f ... in Eq. P^ . 
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FIG. 3: Time dependence of the determinant Ar[47r — 0.6]. The distribution function was taken to be nd(e), Eq. (|27|l . Numerical 
results for matrix sizes N — 50, 100, and 500, as well the leading asymptotic term (|3ip are shown. The numerical result for 
N = 500 is almost indistinguishable from the asymptotics. 
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FIG. 4: Difference between the numerically evaluated deter- 
minant (shown in Fig. [Sjl and the asymptotic approximation 
for N — 5000. The electronic distribution was the double- 
step distribution nd{e). Dotted line: only the leading term 
(|3ip was subtracted; full line: three main harmonics (Fisher- 
Hartwig branches) of the expansion (|15p have been taken into 
account. The remaining difference is due to a correction [of 
the type indicated by -|- . . . in Eq. (fT5|) ] to the leading har- 
monic (1311). 



From now on, we measure t in units of 1/U. Figure [3] 
shows the result of numerical evaluation of the normal- 
ized Toeplitz determinant Ai-[47r — 0.6] for the double- 
stcp distribution function rid given by Eq. (1271) . We have 
plotted the data for N — 50, 100, 500 together with lead- 
ing term of the asymptotic expansion (J15l) , the one with 
uq — ni ~ —1 

AflAir - 0.6] w (0.25 + 0.026i)e-'^^'''+'''''^-^/^^*T^-^'^ . 

(31) 
Here r^ « 77 and eo -I- ei — 1/3. Note the fast conver- 
gence with the increase of the matrix size and perfect 
agreement with the predicted asymptotic behavior. We 
stress that the asymptotic fit used here has no adjustable 
parameters. 

Let us now explore the effect of the other terms in the 



expansion p5p . The next two terms are characterized 
by (no,ni) = (—2,0) and (no,ni) = (0,-2). Apart from 
the exponential damping at scales larger then r^ they de- 
cay as r~°'^^ and r~°-^^ respectively. Since in this case 
powers of the leading and the subleading harmonics are 
substantially different, a reliable observation of the sub- 
leading ones requires more substantial numerical efforts. 
To achieve the required accuracy, we use larger values of 
the matrix size {N = 5000). Note that in Sec. HV^ 
llVCi where we focus on smaller values of the phase shift 
(5, subleading harmonics will be much more pronounced 
and easily seen. 

The difference between the numerically calculated 
Toeplitz determinant (fTU)) and its asymptotic approxi- 
mation [T5]) is shown in Fig. |4l The dotted line cor- 
responds to the difference between the numerical result 
and the leading term (I5T|) . The solid line is the differ- 
ence between the numerical result and the first three 
terms in the expansion (J15p . As expected, inclusion of 
the terms {no,ni) — (—2,0) and (no,ni) — (0,-2) im- 
proves the agreement between the asymptotics and the 
exact results. Indeed, the oscillations at high frequencies, 
that are clearly seen on a dotted line, are absent on the 
solid line. Nevertheless, a clear difference between the 
exact result and the asymptotic formula remains, which 
is predominantly due to the correction [of the type -I- . . . 
in Eq. dTSl) ] to the leading (oscillating with frequency 
ci + (-2 = 1/3) harmonic. 

We have thus demonstrated that, even for a relatively 
large phase (5, the numerical simulations work perfectly 
and that the large-t behavior is fully understood in the 
framework of the asymptotic expansion. In the sequel, 
we will present the results for two physical problems of 
our interest (FES and Luttinger liquid) in the energy 
domain. This is more natural physically (as this corre- 
sponds to spectroscopy measurements) and also gives us 
the possibility to separate the contributions of different 
harmonics in (jl5|) within the same graph. We note that 




FIG. 5: X-ray absorption spectra at different values of the scattering phase So for the double-step distribution function of 
electrons nd(e) specified in Eq. (|27p . The solid lines represent the result of numerical evaluation of Toeplitz determinants 
while the dotted lines show the fits based on the asymptotics (|15II25() . The legend shows the corresponding dephasing time r,^ 
together with the exponents 70 and 71 governing the singular behavior of G^(e) at e = eo = ~l/3 and e — ei — 2/3. 



the Green functions G<{t) are obtained from a Toeplitz 
determinant (or a product of two Toeplitz determinants) 
by multiplication with 1/{At)'^'^^ (with 7 being the zero 
temperature exponent, see the last column of Table |l]). 
Thus, 



G^ie) 



j) G^ie/U); 



(32) 



where the functions G<{e/U) are cutoff independent. 
From now on we omit the energy independent factor 
{U /h.y from the Green functions and measure all the 
energies in units of the characteristic scale U . 



B. Fermi edge singularity 

According to Eq. ([T|), the emission/absorption rates 
out of equilibrium are given by a single Toeplitz determi- 
nant. We analyze the case of a double-step distribution 
function nd(e), Eq. ([27|. first. The results for the differ- 
ent values of the scattering phase 5q are shown in Fig. [S] 

The solid lines in Fig. [5] represent the result of numeri- 
cal evaluation of Toeplitz determinants, while the dotted 
lines show the fits based on the asymptotic formulas (|15p , 
(1251) . Only the dominant terms in the sum (jlSp were re- 
tained (the terms with no = — 1, ni = or vice versa). 
Using the expansion P^ , we are able to calculate the sin- 
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FIG. 6: X-ray absorption spectra of the non-equilibrium FES problem with triple-step distribution of the electrons nt(e) (see. 
Eq. (|28[)'l at relatively small scattering phase So — ±0.2, ±0.4, ±0.8. The solid lines represent the result of numerical evaluation 
of Toeplitz determinants while the dots show the fits based on the asymptotic expansion (|15[) . For the chosen distribution 
function the dephasing rate 1/r^ is identically zero and the singularities are not smeared. The last graph in the second column 
has an inset exemplifying the smooth function SG'' (e) added to the asymptotic expression to fit the numerical data. 
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FIG. 7: X-ray absorption spectra in the non-equilibrium FES problem with triple-step distribution of the electrons. The 
electronic distribution was the same as on Fig. [6] but the phase (5o is now larger. In addition to singularities at Fermi 
edge efc one observes now a singularity at energy eo — ei ± £2 = with the exponent 70-1+2 originating from the term with 
no = —1 , rii — 1 , n2 — —1 in the sum (|15[) . 



gular behavior of G^ (e). The regular part is controlled by 
the behavior of G^ (r) at small r and therefore contains 
the information that is not retained when one uses the 
asymptotic expressions. In order to compare the singular 
behavior predicted by the asymptotic formulas ([T5|) , (^5)) 
with the exact results, we add a smooth function 5G^{e) 



to the Eqs. US]), ([251). We choose (5G>(e) in the form of a 
polynomial of a relatively low order with coefficients that 
are adjusted to optimize the fit. In fact, already a second 
polynomial is sufficient to get a rather good fit, and we 
used it in most of the cases. In several cases we used a 
fourth order poynomial. An example of such a smooth 



10 




FIG. 8: TDOS v-nif-) of left (first column) and right (second column) particles in a non-equilibrium Luttinger liquid. The 
incoming right movers have the double-step distribution nd(e). The incoming left movers are assumed to have zero temperature, 
so that Alt[<5l] = 1. The phase 5r entering the nontrivial determinant A_b[5h] is in the vicnitiy of for the case of VL{e) and 
in the vicinity of 27r in the case of J^is(e), as indicated in the legend. 



function SG-^ (e) is shown in Fig. IH] (see inset of the lower 
right graph). 

In agreement with the analytical predictions, the ab- 
sorption spectra shown in Fig.[5]deinonstrate singular be- 



havior Gppij(e) ~ ( e — fife + 1^ 1 near the Fermi edges 



Cfe , k = 0, 1. Note that the exponents at two edges are 
different, which is a very good demonstration of the im- 
portance of summation over all Fisher-Hartwig branches 
in Eqs. ([T5|) . (|25|) . One observes the enhancement of ab- 
sorption near the Fermi edges the for So > 0. Contrary, 
for (5o < the absorption is suppressed. Upon increase 
of the modulus of the scattering phase Sq, the exponents 
7/j and the inverse dephasing time r7 grow by absolute 
value. Simultaneously, the dephasing increases, which 



induces a stronger smearing of singularities. 

In Fig. m we plot the results for triple-step distribu- 
tion, rit(e), and for relatively small values of the scatter- 
ing phase ^0- At chosen Sq the dominant terms in the 
expansion (TTSl) are those with all rii — except for one 
Uk — —1 and the only visible singularities are located at 
the Fermi edges e^. In contrast to the case of double-step 
distribution, the growth of the scattering phase Sq is not 
accompanied by smearing of the singularities, since the 
dephasing rate l/r^ is identically zero. 

As 6 — 2tt — 2So deviates further from 2tt, additional 
terms in the series (1151) become important, as illustrated 
in Fig. [71 In particular, at Sq = 0.8, apart from the singu- 
larities at Fermi edges efc a new one (at e = eg — ei + £2 = 
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FIG. 9: Non-monotonous dependence of the non-equilibrium dephasing in a Luttinger liquid on the interactions strength. The 
incoming right movers have double-step distribution nd(e). The incoming left movers are assumed to have zero temperature. 
The interaction is now sufficiently strong so that the phase 5r governing the density of states for the left movers is close to 2'k. 
Correspondingly, the dominant singularities in i^L(e) are now located at ei and £2- The exponents 70 and 71 are positive and 
large and the singular behavior of !^L(e) is difficult to see directly (left panels). It becomes evident, however, if one considers 
the second derivative of TDOS, v'lii) (right panels). 



0) with the exponent 70-1+2 is clearly seen. It originates 
from the term uq = — 1 , rii = 1 , n2 = —1 in Eq. (J15l) . 
This once more confirms the extended Fisher-Hartwig 
conjecture (jlSp with the summation over all branches. 



C. Tunneling into non-equilibrium Luttinger liquid 

Let us now turn to another application of Toeplitz de- 
terminants, the tunneling into the Luttinger liquid. We 



begin by considering the simplest case, when the incom- 
ing right-moving electrons have the double-step distri- 
bution function nd(e), while the left movers are held at 
zero temperature. In this case the determinant Airif^L] 
in Eq. (|3]) is identically equal to unity. If the interac- 
tion is not too strong and one is interested in the density 
of states for the left-movers, the phase Sj^ entering the 
non-trivial determinant Aj^ri^n] is close to zero. On the 
other hand, the correlation functions of the right-movers 
are given by the determinants at phase Sr close to 27r. 
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FIG. 10: The density of states in the Luttinger liquid coupled to two reservoirs with the triple-step distribution nt(e). There 
is no dephasing, so that the singularities are sharp. At weak interaction {K = 0.7, upper right plot), the most important 
terms in the right determinant A_rt[5_r] oscillate as a function of time r with frequencies e^, fc = , 1 , 2. On the other hand, 
the leading contribution to the left determinant decays without oscillations, and the oscillatory terms are very small (decay 
with much larger exponents). This leads to the singular behavior of the density of states at ejj. As the Luttinger parameter 
decreases {K = 0.5, upper left plot), the subleading oscillating terms in Alt[5l] come into play. This leads to additional (weak) 
singularities in i^i?(e) at e = —5/4 and e = —1/4. The corresponding regions of energy are magnified in the lower plots. 



Correspondingly, the dominant singularity in the density 
of states vl (e) for the left particles is the one at e = 
while main singularities of i'R{e) are at e = eg , ei. This 
behavior is illustrated in Fig. [8l Note that the left- 
moving electrons are dephased much stronger^ than the 
right-moving. 

The behavior of vl (e) at stronger interaction (see Fig. 
[9|) demonstrates the non-monotonous dependence of the 
dephasing on the Luttinger liquid parameter K. For 
K < {i - V5)/2 « 0.38, the phase 5^ > tt, and the 
leading singularities in v^ are those at eo and ei. They 
can be clearly seen if one plots the second derivative of 
the density of states with respect to energy (Fig. [HI left 
panel). Note that the smearing of those singularities de- 
creases (i.e. singularities sharpen) with increasing inter- 
action strength K = 0.3 — > 0.25 — >■ 0.2, as K evolves 
towards K — 2> — 2-^2 « 0.17, where 5r — 2tt and the 
dephasing is absent. 

Finally, we consider an interacting wire with triple-step 
distribution nt for both left and right moving electrons. 
In this case, both determinants in Eq. ^ are nontriv- 
ial. The corresponding density of states is shown in Fig. 
[TUl At weak interaction (K = 0.7, upper-left panel of 



Fig. ITU)) , the right determinant Aflr[<5ii:] oscillates as 
a function of time r with frequencies efc, k = 0,1,2, 
while the left determinant decays mostly without oscilla- 
tions. This leads to the singular behavior of the density 
of states at efe. As the Luttinger parameter decreases 
{K — 0.5, upper-right panel of Fig fTU]) . sub-leading os- 
cillating terms in Air[(5i] come into play and additional 
(weak) singularities in h'B,{e) appear at e = —5/4 and 
e = -1/4. 



V. SUMMARY AND OUTLOOK 

To summarize, we have explored single-particle Green 
functions were of many-body fermionic systems in non- 
equilibrium settings characterized by multiple-step en- 
ergy distribution functions. By using a periodic ultra- 
violet regularization, the problem is reduced to that of 
Toeplitz determinants. We have carried out numerical 
calculation of the corresponding Toeplitz determinants 
and thus obtained the results for the non-equilibrium 
Green functions in the entire energy range. Further, by 
employing the extended Fisher-Hartwig conjecture, we 



13 



have analytically determined the energy dependence of 
the Green functions near each of the Fermi edges. 

The obtained Green functions show, in the energy rep- 
resentation, power-law singularities near multiple edges. 
The singularities are in general characterized by different 
power-law exponents and are smeared by dephasing pro- 
cesses. In the special case of a distribution function with 
population inversion that alternates between n = 1 and 
n = 0, the dephasing is absent (i.e. the singularities are 
sharp) and the TDOS (or the absorption rate) exhibits 
enhancement and suppression in alternating succession. 

The results of the numerical and analytical methods 
perfectly agree, thus confirming the validity of the ex- 
tended Fisher-Hartwig conjecture. 

We close the paper by listing some of future research 
directions: 

• It would be interesting to see whether an explicit 
form of correction terms within each harmonic 
[those abbreviated by -I- . . . in Eq. ([TS])] can be 
found. Further, a rigorous mathematical proof of 
the extended Fisher-Hartwig conjecture would be 
certainly desirable. 

• One can consider many-body correlation functions 
in the non-equilibrium setups discussed above. 
This problem can be reduced to determinants that 
are of a form more general than the Toeplitz one. 
Some results in this direction will be reported 



SOO] 
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• It would be important to further extend the Fisher- 
Hartwig conjecture in order to include Toeplitz ma- 
trices with matrix symbols. 
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